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ABSTRACT 

The  temperature  dependence  of  diffusion  length  for 
PU-BE  neutrons  in  distilled  water  poisoned  with  boric 
acid  was  studied  using  concentrations  of  1.0,  8.0,  16.0, 
and  30.0  grams  per  liter.   Runs  at  each  concentration  were 
made  at  temperatures  of  95,  131,  and  167°  F.   The  trend  of 
temperature  dependence  was  clearly  established.   A  theo- 
retical temperature  dependence  of  k2,  where  k  =  1/L  is- 
proposed  using  a  Maxwellian  distribution  and  found  to  be 
in  good  agreement  with  the  data  reported  in  Reference  [2]. 
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I.   INTRODUCTION 

The  dependence  of  thermal  neutron  diffusion  length 
upon  poison  concentration  in  water  has  been  the  subject  of 
several  experiments,  the  most  extensive  of  which  was  done 
by  Goddard  and  Johnson  [2].   Using  boron,  a  strong  1/v 
absorber- -  its  cross  section  varies  inversely  with  neutron 
velocity- -  these  investigators  carried  out  experiments  over 
a  wide  range  of  boron  concentrations. 

The  object  of  this  study  was  to  determine  diffusion 
length  as  a  function  of  boron  concentration  and  tempera- 
ture.  An  attempt  was  made  to  represent  both  the  absorp- 
tion and  the  temperature  dependence  of  diffusion  length  in 
a  polynomial  approximation.   This  experiment  was  conducted 
using  a  plutonium-beryllium  source. 


II.   THEORY 


A.    ONE  GROUP  THERMAL  NEUTRON  DIFFUSION 

For  a  point  source  of  thermal  neutrons  in  an  infinite 
hydrogeneous  medium,  the  one  group  thermal  diffusion  equa- 
tion may  be  written  in  spherical  coordinates  as 

i    a  f    d  $T  1    1 

1     Q         2  J_        1    -       n 

r2   dr   I1     dr  I  "  L2  $T   u     (1) 


Usirg  the  source  condition 
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the  solution  to  the  diffusion  equation  is, 


c  -r/L 

3)   =  Se   / 
T    4-nDr 


(2) 


Where  S  is  the  source  strength  in  neutrons/sec. 

The  determination  of  the  diffusion  length  would  now 
be  possible  using  the  above  solution  if  the  neutron  source 
were,  in  fact,  thermal.   However,  the  PU-BE  neutron  source 
emits  neutrons  with  an  average  kinetic  energy  of  approxi- 
mately 4.2  MEV.   The  actual  neutron  spectrum  extends  from 
very  low  energies  to  about  10.6  MEV.   In  order  to  use 
the  one  group  thermal  model,  it  is  necessary  to  account 
for  the  PU-BE  neutron  energy  spectrum.   This  is  accomplished 


by  considering  the  neutrons  as  being  born  in  two  groups ; 
a  fast  group  with  an  energy  range  above  .025  ev  up  to  10.6 
MEV  and  a  thermal  group  with  energy  .0  25  ev .   The  next 
simplifying  assumption  is  to  assume  all  the  neutrons  born 
to  the  fast  group  are  born  at  the  average  group  energy. 
The  actual  energy  distribution  is  shown  in  Figure  1,  hence 
there  would  be  many  different  migration  areas  [6].   It  is 
assumed  that  the  fast  neutrons  are  born  at  an  average 
energy  and  therefore,  there  will  be  an  average  migration 
area.   The  average  energy  of  the  PU-BE  neutrons  is  4.2  MEV. 

Now  having  lumped  the  fast  neutrons  into  a  single 
energy  group,  the  source  term  in  the  diffision  equation  can 
be  derived.   This  term  must  represent  those  source  neutrons 
that  are  born  fast  and  have  slowed  down  to  thermal  energy. 
The  one  group  diffusion  equation  written  with  the  source 
term  has  the  following  form: 

L_  1_   r2  _T  _  1_   &    +  S(r)   =  0        r2A>, 
r2   dT  r   ^F^   L2   *T     D  lZAJ 

B.    SOURCE  TERM  DEVELOPMENT 

For  an  isotropic  point  source  in  an  infinite  medium, 

the  source  releases  S   neutrons  per  unit  time.   Since  in 

o  F 

hydrogeneous  media  a  fast  neutron  may  loose  all  energy 
after  one  collision,  the  density  of  first  scattering 
collisions  is  given  by  [5] . 

q  (r,  ETH)  =  ^  Zs(E0)e-VV  (3) 


For  an  infinite  hydrogen  medium,  the  slowing  down 
density  at  any  energy  is  constant  and  equal  to  the  source 
density  [4] 


S(r)  =  q  (r,E) 


(4) 


Therefore,  for  a  point  source  emitting  fast  neutrons  at 
E  ,  the  number  of  neutrons  entering  the  thermal  region  is 


given  by 


S(r)  =  q  (r,  ETH)  = 
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substituting  into  equation  (2A)  results  in 
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where , 


C  = 


S  E(E  ) 
o  v  oJ 

4ttD 


and 


1/b  =  W 


C.     SOLUTION  OF  DIFFUSION  EQUATION 

In  the  development  a  first  collision  density  is  as- 
sumed and  E  is  assumed  constant  at  the  average  energy  of 
the  fast  neutron.  The  homogeneous  solution  to  this  equa- 
tion is  : 
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E(n) 

FIGURE  1.   Neutron  Spectrum  for  a  PU-BE  Source 
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Since  the  flux  must  be  finite  everywhere,  B  must  be  zero 
Assuming  a  particular  solution  of 

-r/b 


E  e 


YP   ^        r2  (9) 

and  substituting  in  Equation  (7)  yields 


*T 

— 

+    *P 

*T 

= 

A    e 

-r/L 
r 

Cb2L2e 


r/b 


Dr2(L2-b2)         (10) 

This  is  the  asymptotic  solution  (i.e.,  the  solution  for 
r-K»)  for  Equation  (6)  . 

If  L2  is  larger  than  b,  the  first  term  of  Equation  (10) 
becomes  much  more  significant  than  the  second  term  at  some 
sufficiently  large  r  and  the  flux  can  be  assumed  to  decay 
as  a  simple  exponential.   For  PU-BE  sources  b  is  larger 
than  L2 ,  hence  the  complete  expression  must  be  retained 
to  correctly  describe  the  neutron  diffusion.   The  several 
unknows  associated  with  this  equation;  namely,  the  con- 
stant A,  the  parameters  D,b  considerably  increase  the 
difficulty  of  determining  an  accurate  diffusion  length. 
Even  though  it  has  been  possible  to  develop  a  theoretical 
expression  for  flux  based  on  a  one  group  model  with  a 
source  term  to  account  for  the  fast  neutrons,  this  ex- 
pression is  not  simple  and  does  not  lend  itself  to  a 
straightforward  determination  of  diffusion  length.   How- 
ever, it  is  still  possible  to  proceed  by  turning  to  the 
relationship  between  diffusion  length,  age,  and  migration 


11 


area.   Although  this  method  is  far  from  simple  it  does 
offer  a  solution.   The  next  step  will  be  to  develop  the 
expression  for  age  and  migration  area. 

D.    NEUTRON  AGE  CALCULATION 

The  definition  of  age  can  be  derived  from  the  reactor 
critical: ty  equation  written  in  the  moment  form  [6] 


5      (-l)n  B2n  ?2n 

n=0     - *—r~. — rrn =  1 


1+L2B2    n~u         (2n+l) !     '  x  (12) 

For  an  infinite  medium  it  is  possible  to  neglect  all  terms 
greater  than  n  =  1.  The  summation  of  equation  (12)  can  be 
written  as : 

B272     B-T1*    B6?6 
-  +  +  +  •  •  • 

1     3!       5!      7! 

Now  setting  the  above  approximation  equal  to  the  critically 
equation  in  the  form 

K°°  e _  1 

1+  L2B2  (13) 

we  get 

-B2t    ,         B2?2 
e     =  1  ■ 


3! 

when  all  terms  in  the  series  with  n>l  has  been  neglected. 
Using  a  series  expansion  of  the  exponential,  the  expression 
becomes 

1  -  B2t  =  1  - 


3!  (14) 
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which  yields : 


(E   E)  =  r   (Eo,E) 

6  (15) 


The  age  relates  to  the  mean-square  distance  traveled  by  a 
neutron  from  its  point  of  birth  as  a  fast  meutron  to  the 
point  at  which  it  becomes  thermal. 

The  mean-square  (net  vector)  distance,  r2  ,  which  a 
neutron  travels  from  its  source  to  where  it  is  absorbed, 
is  given  by  [3] . 

/  CO 

/   r2   4irr2Z  $  dr 
Y2    =  i a 

o        a  (16) 


Substituting 


rT   E)   =  SL_LL>Ji) 


^    E  Es  (17) 


from  the  age-diffusion  theory,  gives 


-2  .^r  2  q  (r,E)  dV 

r    /    q  (r,E)  dV  (18) 


The  saturated  activity  can  be  written  as 

s       a  (19) 

where  both  the  macroscopic  cross- section  and  the  flux  are 
functions  of  energy.   Integrating  over  the  energy  range  of 
interest  and  substituting  Equation  (17)  into  Equation  (19) 
yields : 
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q  (E)  E  (E)   dE 

A   =  V    — 

s 


5  Es     E       (20) 


For  the  single  resonance  peak  of  indium  it  is  pos- 
sible to  assume:   (1)   the  integral  is  small  except  at  the 
resonance  peak,  (2)  the  slowing-down  density  is  independent 
of  energy  and  (3)  the  moderating  power  (£E  )  is  constant 
over  the  peak.   Equation  (20)  may  now  be  written  as 


As     £Zs     Za^E)   E        (21) 

This  shoi\s  that  the  slowing  down  density  is  directly  pro 
portional  to  the  saturated  activity.  Substituting  Equa- 
tion (2)  into  Equation  (18)  yields: 


—2       r2  A   (r)  dr 
r   =         s  v  J 
s      

As  (r)  dr         (22) 


which  for  cadium- covered  foils  becomes 
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rCD 


. 
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Finally, 
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An  experimental  evaluation  of  Equation  (24)  is  in- 
cluded in  Appendix  B. 

E.    DETERMINATION  OF  MIGRATION  AREA 

The  total  mean-square  distance  a  neutron  travels  from 
its  point  of  birth  to  its  point  of  interaction,  r2  ,  con- 
sists of  two  parts:   (1)   the  mean-square  distance  covered 
in  slowing  down  process,  r2   and  (2)  the  mean-square  dis- 
tance covered  in  the  process  of  thermal  diffusion,  r2  ,  . 
The  first  parameter  correlates  with  the  age,  and  the  second 
parameter  correlates  with  thermal  diffusion  length,  L2  [6]. 
The  migration  area  is  defined  as 

M2  =  t  +  L2  (25) 

The  correlation  between  r  2  and  M2  is  given  by 


ri  -   t 


I 


M" 


oo 

r1*   A  (r)  dr 

2     /  S 


L      ^  Mr 


)  dr 


(26) 


The  integrals  can  be  obtained  in  a  manner  analogous  to  the 
derivation  of  Equation  (24)  for  r  2. 

F.    TEMPERATURE  DEPENDENCE  APPROXIMATIONS 

In  order  to  formulate  the  temperature  dependence  of 

diffusion  length  a  series  expansion  in  absorption  cross 

section  E    of  the  form  suggested  by  reference  [2]  was 
a  o 


utilized. 


K2    =    a0+    cti    £         -a2£2         +a3E3 

u         l      ao  ^      ao  ao 

(27) 
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where  K  is  the  inverse  of  the  diffusion  length,  £   ,  is 

6   '   ao ' 

the  sum  of  the  2200m/sec  macroscopic  absorption  cross  sec- 
tions for  boron  and  hydrogen  and  the  a's  are  coefficients 
to  be  determined.   We  shall  attempt  to  include  the  tempera 
ture  dependence  in  the  expansion  defined  by  Equation  (29) . 

A  first  approximation  of  temperature  dependence  was 
derived  by  assuming  a  Maxwellian  distribution  [4]  for 
neutron  flux. 


*  (E)  =r-^3/2 
nv  *       (ttRT)  ' 


H]    "'   E  e-E/KT 

m 


(28) 


where  n  is  the  density  of  neutrons,  K  is  Boltzmann's  con- 
stant, m  is  the  meutron  mass,  and  T  is  the  temperature  of 
the  medium  in  degrees  Kelvin.   By  setting  the  derivative 
of  Equation  (28)  equal  to  zero  it  can  be  shown  that  the 
most  probably  velocity  is  given  by  [4] 


from  this 


V^  =  1.2  84  Th   x  10. 


Vr 


cm/ 


sec   (29) 


V 


f       \ 

T 

T 

0 

(30) 


where  V   is  the  most  probable  velocity  for  a  Maxwellian 
flux  at  293. 6°K. 

Using  this  relationship  it  is  possible  to  relate  the 
average  thermal  macroscopic  cross  section  as  a  function  of 
temperature  [4] . 
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This  relationship  was  applied  to  the  polynomial  approxi- 
mation with  the  temperature  dependence  being  applied  to  the 
coefficients.   In  order  to  test  the  validity  of  the  approxi- 
mation the  data  reported  by  reference  [2]  at  22.3°C  was  used 
for  the  initial  data  at  T  .   This  approximation  was  found  to 
correlate  very  well  with  the  data  at  65°C  from  the  same 
reference.   The  relatively  large  error  associated  with  the 
present  research  precludes  any  conclusion  more  specific  than 
that  trend  of  the  data  is  in  the  direction  of  the  approxi- 
mation.  Figure  (7)  is  a  plot  of  K2  versus  S   for  various 

a 

temperatures  using  the  temperature  dependent  approximating 
polynomial . 
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III.   DESCRIPTION  OF  THE  EXPERIMENT 

Indium  foil  activation  measurements  were  made  in 
heated  water  poisoned  with  boric  acid  (H-BO,) .   Fifteen 
indium  foils  were  placed  in  a  plexiglass  holder  at  inter- 
vals of  1/2  inch.   The  holder  was  located  in  the  center  of 
a  22  X  30  X  23  inch  tank  made  of  1  inch  plexiglass.   The 
tank  contained  approximately  45  gals,  of  distilled  water. 
-Two,  rheostat  controlled,  500KW  emersion  heaters  were 
used  to  heat  the  tank.   Temperature  control  within  1/2°F. 
was  provided  by  a  Thermistemp  unit  controlling  a  third 
emersion  heater  rated  at  750KW.   The  tank  was  kept  at  a 
uniform  temperature  by  a  variable  speed  electric  stirrer. 

The  tank  was  located  on  top  of  a  second  tank  contain- 
ing the  radioactive  PU-BE  source.   This  steel  tank  was 
2.95  feet  in  diameter,  1.25  high  and  contained  65  gals,  of 
water.   The  sources  were  contained  in  three  stainless  steel 
tubes  arranged  in  a  triangular  configuration..   The  highest 
source  was  located  5/8  inches  from  the  top  of  the  tank. 
The  sources  were  located  directly  under  the  foil  column 
by  three  plexiglass  extension  rods. 

The  three  source  strengths  were  1,  1,  and  3  curies. 
Since  the  sources  were  anisotropic,  the  extension  rods  were 
marked  so  that  the  sources  could  be  placed  to  insure  the 
same  exposure  during  every  run.   Figure  2  shows  the  ex- 
perimental set-up. 


18 


Saturated  boric  acid  solutions  were  added  for  poisoned 
runs.   After  thorough  mixing,  samples  were  drawn  and  an- 
alyzed checmically  as  described  in  Appendix  A. 

Fifteen  foils  were  used  for  each  run.   The  relative 
flux  was  measured  with  bare  indium  foils  1.0  inches  in 
diameter  and  .01  inches  thick.   Indium  resonance  activity 
was  measured  with  similar  foils  placed  in  cadmium  covers 
.020  inches  thick.   Irradiation  time  for  each  run  was  three 
hours.   The  foils  were  removed,  one  at  a  time,  and  allowed 
to  decay  for  three  minutes  before  counting.   This  waiting 
time  allowed  the  13  second  half-life  In     to  decay  to  a 
level  that  would  not  interfere  with  the  54  minute  activity. 
Foils  were  counted  for  three  minutes  using  the  equipment 
specified  in  Table  I.   This  counting  technique  resulted  in 
standard  deviations  ranging  from  .1  to  10  percent.   Back- 
ground counts  were  taken  before  and  after  each  run  and 
averaged.   The  average  background  was  subtracted  from  the 
foil  counts  before  correction  factors  were  applied. 

Detector  GM  tube 

operating  voltage  1300  volts 

Voltage  Supply  Hewlett-Packard  Model  6516A 

Preamplifier  Hewlett-Packard  Model  1045A 

Time-Scalar  Hewlett-Packard  Model  5201L 

TABLE  I. 

COUNTING  EQUIPMENT 
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IV.   ANALYSIS  OF  DATA 

A.    ACTIVATION  DATA 

Each  run  was  conducted  with  fifteen  foils.   The  foils 
were  irradiated  for  three  hours,  removed  one  at  a  time, 
then  counted- starting  with  the  foil  closest  to  the  source. 
In  this  way,  the  foils  with  the  lowest  counts  were  irradia- 
ted as  long  as  possible  to  improve  the  counting  statistics. 
Foil  counts  were  corrected  for  differences  in  foil  weights, 
absorption  of  indium  resonance  neutrons  by  the  cadmium  and 
counter  resolving  time. 

Foil  weights  were  normalized  to  1.0  gram.   Correction 
factors  ranged  from  1.071  to  1.082. 

The  indium  resonance  counting  rate  correction  was 
1.11  as  determined  from  reference  [10], 

Counter  resolving  times  were  determined  by  the  two- 
source  method  using 

Al    =   A2         •    . 


1-N1t     l"N2t 

where  A  is  the  saturated  activity,  N  is  the  mean  counting 
rate,  and  t  is  the  resolving  time.   A  correction  factor 
of  1/(1-Nt)  was  applied  to  all  counting  times. 

The  saturated  activity  was  calculated  using  the  cor- 
rected counts  in  the  equation 

A   -   AC  

(l-e"Atl)  e-^t2"  ^  (1  -  e'^U'U)) 
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where , 

ti  -  irradiation  time 

t2  -  time  counting  started 

t3  -  time  counting  stopped 

t2-ti  -  waiting  time 

t3-t2  -  counting  time 

Having  calculated  the  relative  saturated  activity, 
the  quantities  Asr2  ,  Asr4  ,  A  ,r2  and  A  ^r1*  were  computed 
and  plotted  versus  the  distance  from  the  source.   The  com- 
puter program  used  for  these  calculations  is  described  in 
Appendix  D. 

B.    RELAXATION  LENGTH 

Equation  (6)  gives  the  space  dependence  of  the  slow- 
ing- down  density  as 

C  e-r/b 
q  0,E)  =  — —2 

Since  the  activity  is  proportional  to  the  slowing 

down  density  as  shown  in  Equation  (21) ,  we-  have 

C   e-/b 
A   = 


s       r2  (32) 

Rewriting  and  taking  the  logarithum  of  both  sides  yields 


In  Agr2  =  CONST  -  r/b  .^ 

The  quantity  In  A  r2  was  plotted  against  distance.   The 
relaxation  length  is  the  negative  reciprocal  of  the  slope 
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The  same  method  was  used  to  determine  the  relaxation 
length  associated  with  the  indium  resonance  activity. 

C.  AGE  AND  MIGRATION  AREA 

Equations  (24)  and  (26)  were  used  to  determine  Age 
and  Migration  Area.  Computer  Program  No.  1  (Appendix  D) 
plotted  graphs  of  Asr2  ,  A  j2,  Asr1*  ,  and  A  iTh  versus  r. 
A  smooth  curve  was  drawn  through  the  data  points  up  to  a 
radius  of  12  cm.  At  12  cm  the  plot  was  considered  expo- 
nential and  extrapolated  to  infinity  by  means  of  a  curve 
formed  from  the  equation 

Asr2  =  C  e"r/b 

where  b  is  the  relaxation  length,  determined  from  curve 

fitting  of  measured  activities  and  the  logarithmic  plot. 

The  areas  under  the  Asr2  and  A  ,r2  curves  were  determined 

cd 

by  Simpsons  rule.   The  areas  for  Asr1*  and  A  j11  were  de- 
termined by  multiplying  the  Ar2  data  points  by  r2  and 
integrating  the  result,  rather  than  integrating  the  Asr14  and 
A  ,rk    curves  directly.   After  the  areas  had  been  computed 
the  second  moments  were  determined  by  Equations  (23)  and 
(26) ;  and  finally  Age  and  Migration  Area  were  found  from 
Equations  (24)  and  (26) .   The  computer  program  used  for 
these  computations  is  Appendix  E.   Plots  of  Asr2  and  A  J2 
are  included  in  Appendix  B. 

D.  DIFFUSION  LENGTH 

Diffusion  length  is  determined  from  Equation  (11) . 
The  Age  that  has  been  determined  is  to  the  indium  resonance. 
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1.8  cm2  was  added  to  the  Age  to  account  for  the  Age  in- 
crease from  indium  resonance  to  thermal  energies  [6].   The 
nature  of  the  calculation- the  difference  of  two  relatively 
large  numbers- leads  to  considerable  variation  in  the  re- 
sults even  under  the  best  experimental  conditions.   Small 
errors  (±.5  cm2)  in  the  computed  values  of  the  Migration 
Area  and  Age  lead  to  20  percent  change  in  the  value  of 
diffusion  length.   If  errors  are  estimated  to  be  2.5  cm2, 
a  very  reasonable  estimate  for  this  experiment,  the  value 
of  diffusion  length  at  higher  poison  concentrations  will 
vary  200  percent.   As  the  poison  concentration  is  increased, 
diffusion  length  decreases  and  using  the  same  estimates  of 
error  the  percent  error  becomes  very  large.   From  this  it 
is  concluded  that  the  results  of  this  experiment  serve  only 
to  show  the  trend  of  the  temperature  dependence  of  diffusion 
length  and  the  values  determined  should  not  be  considered 
absolute . 


25 


TABLE  II 


TEMPERATURE  9  5°F 

Boric  Acid  g/1  Migration      Diffusion 

Concentration       Age  Area  Length 


8  ±  .5      58.5  ±  2.1       62.5  ±  1.9      1.5  +  .6 

-  1.3 

16  ±  .6      58.6  ±  1.8       61.6  ±  2.0      1.1  +  .9 

-  1.1 

30  ±  1.1      58.4  ±  1.9       60.9  ±  2.2       .9  +  1.9 

-  .9 


TABLE  III 

TEMPERATURE  131°F 

Boric  Acid  g/1  Migration      Diffusion 

Concentration       Age  Area  Length 


8  ±   .5      59.5  ±  1.8       63.6  ±  1.7      1.5  +   .7 

-  1.5 

16  ±   .6      59.6  ±  1.9       60.9  ±  1.8      1.2  +   .6 

-  1.2 

30  ±  1.1      59.6  ±  2.0       62.2  ±  2.2  .     .9  +  1.1 

-  .9 


TABLE  IV 


TEMPERATURE  16 7° F 


Boric  Acid  g/1  Migration      Diffusion 

Concentration       Age  Area  Length 

8  ±   .5     60.2  ±  1.7       64.4  ±  2.1     1.6  +   .8 

-  1.6 

16  ±   .6     60.4  ±  1.7       63.1  ±  2.0      1.2  +   .6 

-  1.2 

30  ±  1.1      62.9  ±  1.9        60.3  ±  2.2       .9  +  1.0 

-  .9 
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__  SQUARE  DAfA  REF.  [2]  AT  22.3°C 

CURVES  FOR  22.3,  35,  55,  65,  70°C 


74 


.5 


FIGURE  7 
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V.   CONCLUSIONS 

Both  Age  and  Migration  Area  were  found  to  be  relatively 
weak  temperature  dependent  parameters  in  the  temperature 
range  under  investigation.   The  trend  of  the  temperature 
dependence  was  clearly  established  and  in  agreement  with 
the  theory.   Errors  associated  with  the  calculations  of 
these  parameters  were  mainly  due  to  the  statistical  errors 
from  the  weak  neutron  sources  and  the  location  of  the  neu- 
tron origin.   The  products  A  r2  ,  A  r1*  ,  A  ^r2  ,  and  A  .r1* 
were  greatly  influenced  by  the  position  of  the  origin.   The 
triangular  configuration  of  the  sources  made  the  selection 
of  one  origin  for  all  neutrons  difficult.   The  first  origin 
selected  was  the  geometric  center  of  the  three  sources. 
The  validity  of  this  reference  point  was  tested  by  com- 
paring the  relaxation  lengths  obtained  for  pure  water  with 
those  reported  in  references  [6,  9].   The  results  were  not 
in  agreement.   The  best  agreement  was  obtained  by  assuming 
an  origin  on  the  top  surface  of  the  uppermost  (also  the 
strongest)  source.   However,  it  should  be  pointed  out  that 
this  forces  the  assumption  that  all  neutrons  are  born  at 
this  point  when  actually  the  neutrons  from  the  other 
source  are  born  at  a  significantly  different  distance. 

The  diffusion  length  in  this  experiment  was  determined 
by  computing  the  square  root  of  the  difference  between  Age 
and  Migration  Area.  Because  of  the  nature  of  the  calcula- 
tion the  relatively  small  uncertainties  associated  with 
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these  parameters  can  and  did  lead  to  very  large  uncertain- 
ties in  the  diffusion  length.   The  investigator  is  hesi- 
tant to  place  too  much  weight  on  the  final  values  of  dif- 
fusion length.   Rather  than  draw  any  conclusions  as  to 
exact  values  of  diffusion  length  the  results  are  reported 
only  to  show  the  trend  of  the  temperature  dependence  of 
diffusion  length. 

Using  a  polynomial  expansion  suggested  by  Ref  [2]  a 
successful  effort  was  made  in  representing  both  the  temp- 
erature and  the  absorption  dependence  of  diffusion  length. 
The  form  of  the  temperature  dependence  was  derived  from 
the  Maxwellian  distribution  and  applied  to  the  macroscopic 
absorption  cross  section.   The  plots  of  the  approximation 
were  in  excellent  agreement  with  the  data  reported  in 
reference  []].   It  is  hoped  that  this  approach  will  serve 
as  a  basis  for  further  work  in  this  area  and  that  the  ap- 
proach will  prove  successful  at  higher  temperatures  and 
poison  concentration. 
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APPENDIX  A:   DETERMINATION  OF  BORIC  ACID  CONCENTRATION 

A.  1  gram  of  Boric  Acid  (H_BO_)  was  dissolved  in  300  ml 
of  distilled  ELO.  Twenty- five  ml  of  glycerol  and  10  drops 
of  phenophthlen  indicator  were  then  added  and  the  result- 
ing solution  was  titrated  using  a  .IN  NaOH  solution.   A 
distinct  reddish-pink  color  marked  the  titration  end-point. 
This  titration  established  the  molarity  of  the  standard. 

B.  Sample  Analysis 

A  25  ml  sample  was  drawn  for  analysis.   To  this  sam- 
ple were  added  5ml  of  glycerol  and  4  drops  of  phenophthalein 
This  sample  was  then  titrated  using  a  .IN  NaOH  solution. 

C.  Sample  Calculation 

143  ml  of  standard  (.IN)  NaOH  solution  was  used  to 
reach  end-point  of  the  standard  solution  (1  gm  dissolved 
in  300  ml  H-0)  .   27.5  ml  of  NaOH  produced,  end-point  in 
sample. 

2  7.5  ml  y  ml 

25  ml  (SAMPLE)  300  ml  (STANDARD) 

y  =  33  ml 
From  standard 

SOT^T     =    3'33   S^/liter 

143  ml    =  33ml 
3.33  g/1     x 

x   =    7 . 7  gm/liter 
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APPENDIX  B 


LOTS  OF  A  r2  and  A  Ar2    FOR  POISON  RUNS 
s        cd 
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APPENDIX  C 


ACTIVATION  DATA  FOR  POISON  RUNS 
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TABLE  V 

-BO-  CONCENTRATION  8.0 
5      5  TEMPERATURE  9  5°F 

g/1;  INDIUM  RESONANCE  ACTIVATION 

r  (cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

4338.9062 

3510.6030 

59.25 

5.08 

2691.9985 

2172.0639 

46.60 

6.35 

2028.5357 

1631.6850 

40.39 

7.62 

1447. 7858 

1160.9482 

34.07 

8.89 

1094.8874 

8713.8183 

29.51 

10.16 

769.0697 

612.2907 

24.74 

11.43 

580.6441 

460.4277 

21.45 

12.70 

423.2587 

334.1813 

18.28 

13.97 

315.7761 

248.0522 

15.74 

15.24 

235.0387 

183.7874 

13.55 

16.51 

171.5650 

133.4304 

11.55 

17.  78 

127.5265 

98.6032 

9.92 

19.05 

89.6837 

68.9076 

8.30 

20.32 

67.0308 

51.1201 

7.14 

21.59 

75.2695 

40.5087 

6.36 

TABLE  VI 


H  BO-  CONCENTRATION  8.0  g/1;  TOTAL  ACTIVATION 
5     TEMPERATURE  9  5°F 


r(cm)      Ag(CPM) 


3 

5 

6 

7 

8 

10 

11 

12 

13 

15 

16 

17 

19 

20 

21 


81 
08 

35 
62 
89 
16 

43 
70 
97 
24 
51 
78 
05 
32 
59 


11509 

9497 

7423 

5566 

4071 

2972 

2250 

1729 

1235 

989 

732 

522 

419 

322 

248 


8781 
4623 
16  5  7 
3006 
2426 
1105 
4948 
0450 
7078 
6018 
5951 
9899 
0863 
2161 
4556 


MEAN  COUNT 


13952 

8464 

6507 

4894 

3612 

2627 

1981 

1505 

1077 

858 

622 

473 

356 

2  72 

208 


5742 
7929 
8710 
4101 
4257 
4343 
6005 
0419 
4638 
1154 
7670 
8107 
7763 
4084 
4776 


STANDARD 
DEVIATION 


118 
92 
80 
69 
60 
51 
44 
38 
32 
29 


12 
00 
67 
96 
10 
25 
51 
79 
82 
29 


24.95 
21.76 
18.88 
16.50 
14.43 
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TABLE 

VII 

_B0,  CONCENTRATION  8.0 
5           TEMPERATURE  131°F 

g/1;  INDIUM  RESONANCE  ACTIVATION 

r  (  cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

3740.3436 

3040.1044 

55.13 

5.08 

2271.1732 

1841.5061 

42.91 

6.35 

1736.7794 

1404.8835 

37.48 

7.62 

1290.1232 

1041.6748 

32.25 

8.89 

889.0577 

715.0063 

26.73 

10.16 

675.6499 

541.6247 

23.27 

11.43 

540.5715 

431.8361 

20.78 

12.70 

365.6290 

290.9912 

17.05 

13.97 

269.6938 

213.8565 

14.62 

15.24 

207.0582 

163.2117 

12.77 

16.51 

159.8396 

124.7460 

11.16 

17.78 

110.9198 

86.0412 

9.27 

19.05 

107.6441 

83.0443 

9.11 

20.32 

72.8171 

55.8133 

7.47 

21.59 

60.8034 

46.2799 

6.80 

11 


TABLE 

VIII 

-BO_  CO?> 

5        5  ra 

ICENTRATION  8.0 

g/1;  TOTAL  ACTIVATION 

IPERATURE  131°F 

r(cm) 

AS(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

20109.0553 

18046.6484 

134.33 

5.08 

11566.5587 

10350.9726 

101.73 

6.35 

8869.6275 

7913.3593 

88.95 

7.62 

6536.7874 

5812.9609 

76.24 

8.89 

4962.5813 

4397.5390 

66.31 

10.16 

3720.1893 

3284.1132 

57.30 

11.43 

2812.3892 

2472.5913 

49.72 

12.70 

2091.6808 

1830.8740 

42.78 

13.97 

1610.7691 

1403.2451 

37.45 

15.24 

1271.6344 

1101.6093 

33.19 

16.51 

935.2484 

806.1274 

28.39 

17.78 

729.9680 

625.0952 

25.00 

19.05 

546.4106 

465.1699 

21.56 

20.32 

412.6101 

348.8295 

18.67 

21.59 

327.0331 

274.4113 

16.56 

55 


II 


TABLE 

IX 

~B0_  CONCENTRATION  8.0  g 
5      *   TEMPERATURE  16 7° F 

/l;  INDIUM  RESONANCE  ACTIVATION 

r  (cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

4268.5135 

3436.8632 

58.62 

5.08 

2863.6486 

2299.7617 

47.95 

6.35 

2335.1988 

1870.2941 

43.24 

7.62 

1668.7725 

1332.7524 

36.50 

8.89 

1369.9945 

1090.5432 

33.02 

10.16 

973.7141 

772.9123 

27.80 

11.43 

697.1962 

551.6057 

23.48 

12.70 

475.3571 

374.7990 

19.35 

13.9  7 

388.4696  ' 

303.0793. 

17.40 

15.24 

392.8802 

255  .2706 

15.97 

16.51 

264.9628 

205.1085 

14.32 

17.76 

316.3616 

162.4721 

12.74 

19.05 

139.1435 

106.7854 

10.33 

20.32 

107.4450 

82.6591 

9.09 

21.59 

89.9828 

68.8904 

8.30 

II 


TABLE 

X 

-BO.,  CONCENTRATION  8.0  g 
6          TEMPERATURE  16 7° F 

/l;  TOTAL  ACTIVATION 

A        /*  f~\    T"\  A  *   N 

STANDARD 

r  (cm) 

Ag(CPM) 

MEAN  COUNT 

DEVIATION 

3.81 

20692.8863 

18570.5976 

136.27 

5.08 

12481.0852 

11170.0156' 

105.68 

6.35 

9169.2326 

8180.6640 

90.44 

7.62 

6277.5951 

5582.4726 

74.71 

8.89 

5344.6326 

4736.0898 

68.81 

10.16 

3937.0296 

3475.5358 

58.95 

11.43 

2886.7369 

2537.9560 

50.37 

12.70 

2290.5178 

2004.9187 

44.77 

13.97 

1768.2372 

1540.4257 

39.24 

15.24 

1368.2429 

1185.8757 

34.43 

16.51 

1108.0042 

955.0327 

30.90 

17.78 

861.6697 

738.2927 

27.17 

19.05 

659.6030 

561.5327 

23.69 

20.32 

477.2240 

403.4556 

20.08 

21.59 

346.5687 

290.9023 

17.05 

56 


TABLE  XI 


H  BO   CONCENTRATION  16  g/1;  INDIUM  RESONANCE  ACTIVATION 
6      *   TEMPERATURE  9  5°F 


r(cm) 


ACD(CPM) 


MEAN  COUNT 


STANDARD 
DEVIATION 


3 

5 

6 

7 

8 

10 

11 

12 

13 

15 

16 

17 

19 

20 

21 


81 
08 

35 
62 
89 
16 
43 
70 
97 
24 
51 
78 
05 
32 
59 


4116 

3034 

2185 

1419 

1164 

755 

558 

455 

309 

227 

220 

209 

131 

130 

117 


3652 

3682 
4791 
0948 
1740 
370  4 
3584 
3512 
0079 
6  5  72 
7544 
2343 
9818 
0777 
1275 


3151 

2334 

1689 

1101 

90  7 

591 

438 

359 

244 

180 

175 

167 

105 

104 

94 


4687 
3879 
1000 
6096 
5075 
1831 
6540 
0251 
4800 
7092 
7788 
1031 
7061 
4638 
3071 


56 
48 
41 
33 
30 
24 
20 
18 
15 
13 
13 
12 
10 
10 
9 


13 

31 
09 
19 
12 
31 
94 
94 
63 
44 
25 
92 
28 
22 
71 


H 


TABLE 

XII 

BO^  CONCENTRATION  16  g/1 
5   4  TEMPERATURE  9  5°F 

;  TOTAL  ACTIVATION 

r  (cm) 

Ag(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

18626.1055 

9497.1953 

97.45 

5.08 

10895.5391 

5582.4882 

74.71 

6.35 

6685.8054 

3441.4116 

58.66 

7.62 

5595.1975 

2890.9396 

53.76 

8.89 

3702.6471 

1922.2897 

43.84 

10.16 

2939.2058 

1532.0256 

39.14 

11.43 

2357.4266 

1233.4499 

35.12 

12.70 

9306.7335 

804.2429 

28.35 

13.97 

1391.1851 

733.0488 

27.07 

15.24 

970.9109 

513.2778 

22.65 

16.51 

728.9893 

386.5915 

19.66 

17.78 

516.2110 

274.5700 

16.57 

19.05 

404.9721 

216.0157 

14.69 

20.32 

313.6666 

167.7666 

12.95 

21.59 

230.3158 

123.5047 

11.11 

57 


II 


TABLE 

XIII 

_B0^  CONCENTRATION  16  g/ 
5   a  TEMPERATURE  131°F 

1;  INDIUM  RESONANCE 

ACTIVATION 

r(cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

4133.0768 

3164.2629 

56.25 

5.08 

2710.5077 

2085.2368 

45.66 

6.35 

2175.8190 

1681.6340 

41.00 

7.62 

1594.8277 

1238.0266 

35.18 

8.89 

1188.0543 

926.1230 

30.43 

10.16 

876.8356 

686.2463 

26.19 

11.43 

687.5132 

540.1201 

23.24 

12.70 

501.3977 

395.3308 

19.88 

13.97 

382.6678 

302.7580 

17.39 

15.24 

282.6519 

2  2  4.3627 

14.97 

16.51 

224.8465 

179.0  371 

13.38 

17.78 

181.8621 

145.2425 

12.05 

19.05 

127.9392 

102.4683 

10.12 

20.32 

97.5775 

78.3633 

8.85 

21.59 

73.5245 

59.1995 

7.69 

TABLE  XIV 


H-BO-  CONCENTRATION  16  g/1;  TOTAL  ACTIVATION 
3  TEMPERATURE  131°F 


r(cm) 

Ag(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

15732.0113 

8021.5351 

89.56 

5.08 

9167.6452 

4697.1757 

68.53 

6.35 

6695.2128 

3446.2539 

5  8.70 

7.62 

5151.0925 

2663.1140 

51.60 

8.89 

3788.3850 

1966.8022 

44.34 

10.16 

2832.8391 

1476.5834 

38.42 

11.43 

2144.8511 

1122.2265 

33.49 

12.70 

1595.8782 

838.0163 

28.94 

13.97 

1177.5317 

620.4697 

24.90 

15.24 

906.2015 

479.0686 

21.88 

16.51 

668.0390 

354.2690 

18.82 

17.78 

505.9825 

269.1296 

16.40 

19.05 

354.5358 

189.1125 

13.75 

20.32 

287.2258 

155.6246 

12.39 

21.59 

213.6016 

114.5419 

10.70 
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TABLE 

XV 

H..BO.,  CONCENTRATION  16  g/1 
5           TEMPERATURE  16 7° F 

;  INDIUM  RESONANCE 

ACT I VAT IO 

STANDARD 

r  (cm) 

ACD(CPM) 

MEAN  COUNT 

DEVIATION 

3.81 

4213.6345 

3225.9375 

56.79 

5.08 

2672.1145 

2055.7004 

45.33 

6.35 

2102.6545 

1625.0871 

40.31 

7.62 

1520.5024 

1180.3295 

34.35 

8.89 

1145.9158 

893.2749 

29.88 

10.16 

869.8834 

680.8051 

26.09 

11.43 

673.7225 

529.2858 

23.00 

12.70 

492.3515 

388.1984 

19.70 

13.97 

365.0609 

288.8278 

16.90 

15.24 

278.6658 

221.1986- 

14.87 

16.51 

211.2086 

168.1763 

12.96 

17.78 

158.9448 

126.9  39  8 

11.26 

19.05 

123.8967 

99.2307 

9.96 

20.32 

90.8913 

72.9937 

8.54 

21.59 

70.3152 

56.6154 

7.52 

TABLE  XVI 


H-BO.,  CONCENTRATION  16  g/1;  TOTAL  ACTIVATION 
*      *   TEMPERATURE  16  7°F 


r( 

;cm) 

Ag(CPM) 

3, 

,81 

15683.6585 

5. 

,08 

9559.2448 

6. 

,35 

7208.4950 

7, 

,62 

5502.6384 

8. 

,89 

3895.1716 

10. 

,16 

2960.5430 

11, 

,43 

2208.6044 

12, 

,70 

1682.5217 

13, 

,97 

1226.8919 

15, 

,24 

936.2943 

16, 

,51 

698.5095 

17, 

,78 

500.3817 

19, 

,05 

395.2241 

20, 

,32 

299.4661 

21, 

.59 

216.0595 

MEAN  COUNT 


STANDARD 
DEVIATION 


7996 

4897 

3710 

2844 

2022 

1543 

1155 

883 

646 

494 

370 

266 

210 

160 

115 


8828 
8164 
4577 
8632 
2421 
1474 
5834 
5139 
4787 
9772 
4277 
1506 
8160 
1714 
8599 


89 
69 
60 
53 

44 
39 


42 
98 
91 
33 
96 
28 


33.99 
29.72 


25 
22 
19 
16 
14 
12 
10 


42 
24 
24 
31 
51 
65 
76 
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TABLE  XVII 


H_BO^  CONCENTRATION   30  g/1;  INDIUM  RESONANCE  ACTIVATION 
5      5  TEMPERATURE  9  5°F 


r  (cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

3814.3995 

1752.1708 

41.85 

5.08 

2444.9645 

1128.5703 

33.59 

6.35 

1898.0324 

880.1640 

29.66 

7.62 

1466.2170 

682.9135 

26.13 

8.89 

1152.9376 

539.2492 

23.22 

10.16 

807.3318 

379.1101 

19.47 

11.43 

595.1468 

280.5332 

16.74 

12.70 

455.3512 

215.4151 

14.67 

13.97 

351.5190 

166.8683 

12.91 

15.24 

238.8134 

113.7388 

10.66 

16.51 

181.2051 

86.5722 

9.30 

17.78 

125.2512 

60.0184 

7.74 

19.05 

88.8686 

42.7057 

6.5  3 

20.32 

44.0997 

21.2496 

4.60 

21.59 

17.3741 

8.3934 

2.89 

TABLE 

XVIII 

H_BO^  CONCENTRATION   30  g/1;  TOTAL  ACTIVATION 
5  TEMPERATURE  95° F 

r  (cm) 


3 

5 

6 

7 

8 

10 

11 

12 

13 

15 

16 

17 

19 

20 

21 


81 
08 

35" 
62 
89 
16 
43 
70 
97 
24 
51 
78 
05 
32 
59 


Ag(CPM) 


9512 

5909 

4424 

3384 

2552 

1875 

1474 

1058 

773 

579 

422 

322 

242 

202 

110 


9576 
7507 
9  6  2  5 
5104 
1833 
4493 
2346 
7  32  5 
8019 
9679 
9913 
3804 
7463 
5797 
3097 


MEAN  COUNT 


4850 

3027 

2277 

1749 

1525 

965 

767 

556 

408 

307 

225 

171 

129 

108 

59 


5273 
9460 
6787 
7912 
0078 
1225 
2319 
7866 
3164 
0202 
1728 
6843 
6355 
4727 
2159 


STANDARD 

DEVIATION 

69, 

.64 

55, 

,02 

47, 

.72 

41, 

,83 

36, 

.40 

31, 

,06 

27, 

,69 

23, 

,59 

20, 

,20 

17, 

,52 

15, 

,00 

13, 

.10 

11, 

.38 

10 

.41 

7, 

,69 
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II 


TABLE 

XIX 

_B0_  CONCENTRATION  30  g/1 
5      *   TEMPERATURE  131°F 

. ;   INDIUM  RESONANCE  ACTIVATION 

r  (cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

4032.5113 

3048.0051 

55.20 

5.08 

2261.6839 

2021.6328 

44.96 

6.35 

2034.2463 

1574.2250 

39.6  7 

7.62 

1441.6142 

1120.4519 

33.47 

8.89 

1119.1448 

871.3825 

29.51 

10.16 

845.8544 

660.1320 

25.69 

11.43 

639.5929 

501.1230 

22.38 

12.70 

483.9967 

380.6340 

19.50 

13.97 

355.0934 

280.2563 

16.74 

15.24 

272.6566 

215.9255 

14.69 

16.51 

217.1437 

172.5205 

13.13 

17.78 

161.9829 

129.0929 

11.36 

19.05 

120.0964 

95.9932 

9.79 

20.32 

93.7458 

75.1415 

8.66 

21.59 

73.6595 

59.1995 

7.69 

TABLE 

XX 

H  BO.,  CONCENTRATION  30  g/1;  TOTAL  ACTIVATION 
5  *    TEMPERATURE  131°F 


r(cm) 

AS(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

9348.0959 

4766.4648 

69.03 

5.08 

5167.0679 

2647.4218 

51.45 

6.35 

4059.7207 

2089.6765 

45.71 

7.62 

3025.5249 

1554.2170 

39.42 

8.89 

2255.1747 

1171.4912 

34.22 

10.16 

1693.3000 

883.1005 

29.71 

11.43 

1280.2822 

670.5717 

25.89 

12.70 

969.5395 

509.6274 

22.51 

13.97 

733.6620 

386.9528 

19.67 

15.24 

572.7375 

303.0563 

17.40 

16.51 

431.5252 

229.0411 

15.13 

17.78 

323.7266 

172.3308 

13.12 

19.05 

242.8408 

129.6355 

11.38 

20.32 

187.9917 

100.6241 

10.03 

21.59 

138.5682 

74.3591 

8.62 
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TABLE  XXI 


H-BO-  CONCENTRATION  30  g/1;  INDIUM  RESONANCE  ACTIVATION 
6      5  TEMPERATURE  16 7° F 


H 


r(cm) 

ACD(CPM) 

MEAN  COUNT 

STANDARD 
DEVIATION 

3.81 

3956.9440 

3029.4165 

55.04 

5.08 

2678.0214 

2060.2446 

45.38 

6.35 

2029.5357 

1568.5754 

39.60 

7.62 

1547.6544 

1201.4069 

34.66 

8.89 

1174.0066 

915.1723 

30.25 

10.16 

889.3505 

696.0410 

26.38 

11.43 

684.7551 

537.9533 

23.19 

12.70 

519.7670  • 

409.8142 

20.24 

13.97 

393.3046 

310.2597' 

17.61 

15.24 

306.2453 

242.2943 

15.56 

16.51 

239.2349 

189  .8993 

13.78 

17.78 

175.6441 

139.8590 

11.82 

19.05 

137.7631 

110.0232 

10.48 

20.32 

107.9603 

85.8811 

9.26 

21.59 

80.4186 

64.5832 

8.03 

TABLE 

XXII 

^BO^  CONCENTRATION  30  g 
5   5  TEMPERATURE  16 7° F 

/l;   TOTAL  ACTIVATION 

STANDARD 

r  (cm) 

AS(CPM) 

MEAN  COUNT 

DEVIATION 

3.81 

8459.7042 

4365.5546' 

66.07 

5.08 

5483.4342 

2841.7880 

53.30 

6.35 

4247.0050 

2209.9780 

47.01 

7.62 

3191.6215 

1656.6054 

40.70 

8.89 

2346.1385 

1230.7382 

35.08 

10.16 

1762.6128 

927.8654 

30.46 

11.43 

1344.7866 

710.2744 

26.65 

12.70 

997.1895 

528.3569 

22.98 

13.97 

742.7238 

394.7204 

19.86 

15.24 

588.4379 

313.6276 

17.70 

16.51 

440.6199 

235.4891 

15.34 

17.78 

91.4371 

48.9967 

6.99 

19.05 

242.5287 

130.2845 

11.41 

20.32 

200.2300 

107.8186 

10.38 

21.59 

141.4392 

76.3347 

8.73 
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APPENDIX   D 


COMPUTER  PROGRAM  NO.  1 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  (1)  CORRECT  THE 
COUNTING  DATA  FOR  NORMALIZED  FOIL  WEIGHT,  ABSORPTION 
OF  INDIUM  RESONANCE  NEUTRONS  BY  THE  CADMIUM  COVER  AND 
COUNTER  RESOLVING  TIME,  (2)  COMPUTE  THE  MEAN  COUNTING 
RATE  AND,  FROM  THIS,  DETERMINE  THE  STANDARD  DEVIATION, 
(3)  COMPUTE  THE  PRODUCTS  AS*R**2 , AS*R**4, LOG  AS*R**2, 
AS(CD)*R**2,  AS(CD)*R**4,  AND  (4)  PLOT  GRAPHS  OF  THE 
ABOVE  PRODUCTS  VERSUS  R,  THE  DISTANCE  FROM  THE  SOURCE. 


PROGRAM  NOTATION 


Dl  DECAY  CONSTANT  FOR  1NDIUM-116 

T  COUNTER  RESOLVING  TIME 

POS  DISTANCE  FROM  SOURCE 

X  FOIL  IRRADIATION  TIME 

Y  WAI TING  TIME 

Z  WAITING  TIME  PLUS  COUNTING  TIME 

W  CORRECTION  FACTOR  FOR  FOIL  WEIGHT 

CI  FOIL  COUNTS 

SAT  SATURATED  ACTIVITY 

SIGMA  STANDARD  DEVIATION 

AMEAN  MEAN  COUNTING  RATE 

F  CADMIUM  CORRECTION  FACTOR 


IMPLICIT  REAL*8( A-H,0-Z) 

WRITE (6,1) 

FORMAT (T5, 'DISTANCE' ,T25, ' SAT  ACTIVITY* ,T5  6, 'MEAN' , 
1T87, 'SIGMA'//) 

REAL  AC,AD,AA,BB,DISTA,DISTC,B,SUM,C,D,E?U,V,D1»C1,B1 

DIMENSION  ACf25),AD(25),AA(25),BB(25),DISTA(25),DISTC( 

25),  B(3),SUM 
1(4),C(3,3),0(9),E(3) ,U(402 )  ,V(401) 

REAL  LABEL1/4H     /,LABEL2/4H     /,LABEL3/4H     /,LABE 

L4/4H     / 

REAL* 8  I  TITLE (12) 

RUN  NUMBER  THIRTEEN 

COMPUTATION  OF  SATURATED  ACTIVITY  OF  INDIUM  FOILS  AND 
PRODUCTS  AS*R**2  AND  AS*R**4 

Dl=1.28D-2 
T=6.GD-6 
D050  N=l»15 

READ(5,2)  POS,Cl,X,Y,Z 
FORMAT (5F9.0) 
READ (5, 6)  W 
i  F0RMAT(1F9.0) 
DIST=(P0S-1.750)*0.2  54D01 
Cl=W*3.0*Cl/< 1.0DO-Cl*T) 
ASAT=Dl*Cl/( ( 1.0DO-DEXP(-Dl*X) ) *( DEXP( -D1*Y )-DEXP (-D1* 
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16,6X,D25. 16) 


Z)  )  ) 

SAT=ASAT 

AMEAN=Cl/3.0 

SIGMA=AMEAN**0.5D0 

F=1.11D0 

WRITE^6^1  DI ST, FSAT,AMEAN, SIGMA 
4  F0RMAT(T6,F5.2,6X,D25.16,6X,D25. 
A=FSAT*DIST**2.0D0 

WRITE  (6t8)  DIST,A,B1 

DISTC(N)=DIST 

AC(N)=A 

A0(N)=B1 

50  CONTINUE 

rnuDiiTATinN  HF  SATURATED  ACTIVITY  OF  CADMIUM  COVERED 
FOILS  AND  THE  PRODUCTS  AS(CD1*R**2  AND  AS!CD)*R**4 

D060  N=l,15 

READ(5,5)  P0S,C1,X,Y,Z 
5  FORMAT ( 5 F9.0) 

RE AD (5, 7)  W 
7  FORMAT! 1F9.0)        „c,™, 

DIST=(P0S-1. 750  1*0.254001 

SS?SiSci/(li:888=SEXPC-Dl*X>l*(DEXPC-Dl*Y»-DEXPC-Dl* 

Z)  )  ) 

SAT=ASAT 
AMEAN=Cl/5.0 

SIGMA=AMEAN**0.5D0   ,„_...  ct~MA 
WRITF(6,3)  DIST, SAT, AMEAN, SIGMA 
3  F0RMAT(T6,Fb.2,6X,D2  5.16,6X,D2  5.16,6X,D2b.l6) 

A=SAT*DIST**2.0D0 
B1=SAT*DIST**4.0D0 

WRITE (6,9)  DIST,A,B1 
9  FORMAT!/, T12,F5. 2, 6X, 
DISTA(N)=DIST 
AA(N)=A 
BB(N)=B1 

PLOTTING  OF  COMPUTED  PRODUCTS  VERSUS  R 

60  CONTINUE  , , t  _  .  , n . 

READ  (5,30)  (ITITLE(I)  ,1=1,12) 

30  [^rADRAW8(15,DISTC,AC,0,3,LABEL2,ITITLE,0,0,0,0,0,0,6 

REAo'lsfSl]  (ITITLE(I),I  =  1,12)  .   ' 

31  ^LLADRAW8(15,DISTC,AD,0,3,LABEL2,ITITLE,0,0,0,0,0,0,6 

REA0f(5?321     (ITITLE!I),I=1,12) 

32  ^ALLADRAW8(15,DISTA,AA,0,4,LABEL3,ITITLE,0,0,0,0,0,0,6 

READT(5^33)  (ITITLE(I),I=1,12) 

33  CALL ADRAW8(15,DISTA,BB,0, 5, LABEL4, I  TITLE, 0,0,0, 0,0, 0,6 
.9,1, LAST) 


D25.16,6X,D25.16) 


COMPUTATION/PLOTTING  OF  LOG  AS*R**2  AND  LOG  AS(Cg)R**2 
VERSUS  R  AND  COMPUTATION  OF  A  LEAST  SQUARE  hi  I  ur 

wn?F^EFnRATH^  PURPOSE  OF  LEAST  SQUARE  FIT  ONLY  DATA 
JOINTS  FOUR  THROUGH  FIFTEEN  ARE  USED  SINCE  THIS  IS 
THE  ONLY  PART  OF  THE  CURVE  USED  TO  ESTIMATE  THE 


RELAXATION  LENGTH 


D070  N=l, 15 
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AC(N)=ALOG(AC(N) ) 
AA(N)=ALOG(AA(N) ) 
70  CONTINUE 

READ  (5,34)  (  IT  ITLE (  I  )  ,  1=1 ,  12 ) 

34  F0RMAT(6A8) 
B(  1J  =  0.0 
3(2  )  =0.0 
B(3)=0.0 
D025  1=4,15 
B(1)=B(1HAC( I ) 

B(2)  =  B(2H-AC(I)*DISTC(  I) 
B(3»=B(3)+AC( I J*OISTC(  I)** 2 

25  CONTINUE 
D026  L=l,4 
SUM (L) =0.0 

26  CONTINUE 

DO  35  L=l,4 
DO  35  1=4, 15 
SUM(L)=SUM(L)+DISTC( I ) ** ( L  ) 

35  CONTINUE 
C(l,l)=12 
D045  M=l,3 
D045  N=l,3 

IF(M+N.EQ.2)  GO  TO  45 
C(M,N)=SUM(M+N-2 ) 

45  CONTINUE 

WRITE(b,55)  B(l ) ,B(2),B(3) 

55  F0RMAT(2X, 'VALUES  OF  B ' , 10X , 3F 1 5.4 ) 
WR ITE ( 6 , 65 ) ( SUM ( L ) , L  =  l , 4  J 

65  FORMAT  (2X,  'SUM  OF  X  I  •  ,  1.0X  ,  4F  15  .  4  ) 
WR I T L ( 6 , 7  5 ) ( ( C (  M  ,  N  )  , M= 1 , 3 )  , N=l , 3 ) 

75  F0RMAT(2X, 'VALUES  OF  C ■ , 5X , 9F 12 . 4) 
CALL  SIMQ(C,B,3,KS) 
WRITE(6,35)  B(  1)  ,B(2)  ,B(3) 

85  FORMAT (2X, 'SOLUTION' ,6X,3F10.4) 

U  (  1  )  =  1  0  .  0 

DO  9  5  1=1,401 

V( I)=B(1)+U( l)*(B(2)+U(I )*B(3)  ) 

U( 1+1) =U( I J+0.05 
95  CONTINUE 

CALL  DRAW  ( 15 , DI STC , AC , 1 , 3 , LABEL2 , I T ITLE , 0 , 0 , 0, 0, 0 , 0 , 9 

,15,1, LAST) 

CALL  DRAW(401,U,V,2,0,LABEL2,ITITLE,0,0,0,0,0,0,6,9,1, 

LAST) 

B(1)=0.0 

B(2)=0.0 

B(3J=0.0 

D023  1=4,15 

B( 1)=B(1)+AA( I ) 

B(2)=B(2)+AA( I)*DISTA( II 

B(3)=B(3)+AA( I )*DISTA( I)** 2 
23  CONTINUE 

D027  L=l,4 

SUM(L)=0.0 

27  CONTINUE 
D03b  L=l,4 
D036  1=4,15 
SUM(L)=SUM(L)+DISTA( I  I ** { L  ) 

36  CONTINUE 
C(l,  1)  =  12 
D046  M=l,3 
D046  N=l,3 

IF(M+N.EQ.2)  GO  TO  46 
C(M,N) =SUM(M+N-2) 

46  CONTINUE 

WRITE(6,56)  B( 1) ,B(2) ,B(3) 

56  F0RMAT(2X,«VZLUES  OF  B  '  , 10X , 3F 1 5 .4 ) 
WRITE (6, 66) ( SUM(L) ,L=1 ,4) 

66  FORMAT (2X, 'SUM  OF  X  I « , 10X, 4F 1 5. 4 ) 
WRITE ( 6,76) ( (C(M,N) ,M=1,3) ,N=1,3) 

76  FCRMAT(2X, 'VALUES  OF  C • , 5X , 9F 12 . 4) 
CALL  SIMQCC, 3,3,KS ) 
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WRITE(6T86)  B( 1) ,B(2) ,B(3) 
86  FORMAT (2X,  • SOLUT I  ON '  , 6X, 3F 1 0. 4 ) 

U(l  )  =  10.0 

D096  1=1,401 

V(  I  )=B(  1)+U( I )*(B( 2)+U(  I  )*B(3)  ) 

U( 1+1 )=U(I J+0.05 
96  CONTINUE 

CALL  DRAW(401,U,V,2,0,LABEL3,ITITLEV0,0,0,0T0,0,6,9,1, 

LAST) 

CALL  DRAW  ( 1 5, D I STA , AA  ,3 ,4 , LABEL3, I Tl TLE » 0 ,0 , 0, 0, 0 , 0 , 6 

,9, 1,LAST) 

STOP 

END 
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APPENDIX  E 


COMPUTER  PROGRAM  NO.  2 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  COMPUTE  DIFFUSION 
LENGTH.   THE  AREAS  UNDER  THE  CURVED  PORTIONS  OF  THE 
CURVES  AS*R**2,AS*R**4,AS(CD)*R**2  AND  AS(CD)*R**4  ARE 

COMPUTED  BY  SIMPSON'S  RULE.   THE  AREAS  UNDER  THE 
LINEAR  PORTIONS  OF  THE  CURVES  ARE  COMPUTED  BY  THE 
EQUATION  Y=C2*EXP(-R/B).   THE  TOTAL  AREAS  ARE  USED  TO 
SOLVE  FOR  MIGRATION  AREA  AND  AGE  USING  EQUATIONS  (26) 
AND  (24).   DIFFUSION  LENTH  IS  SOLVED  BY 

DL=(MIG-(AGE+1.8)>**.5  WHERE  THE  L.8  CORRECTION  FACTOR 
CORRECTS  FOR  ADDITIONAL  AGE  REQUIRED  TO  REACH  THERMAL 
ENERGIES. 


PROGRAM  NOTATION 


B     RELAXATION  LENGTH 

BCD   RELAXATION  LENGTH  -  CADMIUM  COVERED  FOILS 
R     DISTANCE  FROM  SOURCE 
S2    AREA  UNDER  AS*R**2  CURVE 
SCD2  AREA  UNDER  AS(CD)*R**2  CURVE 
S4    AREA  UNDER  AS*R**4  CURVE 
SCD4  AREA  UNDER  AS(CD)*R**4  CURVE 
MIG   MIGRATION  AREA 
AGE   AGE 

DL    DIFFUSION  LENGTH 

REALMS  H,SUM2,SUM4,SCD2,S4,  SCD4, A2 , ACD2 , AS2 , AS4, ASCD2 
1S2,Y, YCD,C2,R 
DIMENSION  Y(IOO) ,YCD( 100),R(100) 

RUN  NUMBER  SIX  -  CONC  8  GM/L  -  TEMP  131F 

READ(5,10)  ( Y(  I)  ,1=1,10) 

10  FORMAT(IFIO.O) 

B=11.4 

C2=3000000.0 

READ<5,50)  (YCD( I) ,1=1,10) 
50  FORMAT(IFIO.O) 

BCD=9.0 

CD2=710000.0 

COMPUTING  SATURATED  ACTIVITY  TIMES  DISTANCE  SQUARED  AR 
AS*R**2 

D0160  1=11,80 

R(  I )=1.0*( I) 

Y(  I  )=C2*DEXP(-R( I >/B) 

H=1.0 

N=79 

SUM2=0.0D0 
SUM4=0.0D0 
M=N~1 
160  CONTINUE 

D020  1=3, N, 2 
SUM2=SUM2+Y( I ) 
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20  CONTINUE 

D030  1=2, M, 2 

SUM4=SUM4+Y(  I  ) 
30  CONTINUE 

S2=H/.3D1*(Y(I )+.401*SUM4+.2Dl*SUM2+Y(N+l) ) 

WRITE(6T40)  S2 
40  F0RMAT(T19, 'SIMPSONS  APPROXI MATI ON' , D25. 16/// ) 

AS(CD)*R**2 

D0230  1=11,80 

R(i)=i.o*m 

YCD( I )=CD2*DEXP(-R(I )/BCD) 

H=1.0 

N=79 

SUM2=0.0D0 

SUM4=0.0D0 

M=N-1 
230  CONTINUE 

D060  1=3, N, 2 

SUM2=SUM2+YCD( I ) 
60  CONTINUE 

D070  1=2, M, 2 

SUM4=SUM4+YCD( I ) 
70  CONTINUE 

SCD2=H/.3D1*(YCD( I )+.4Dl*SUM4+. 2D1*SUM2+YCD ( N+l ) ) 

WRITE(6,80)  SCD2 
80  FQRMAT(T19,' SIMPSONS  APPROXIMATION' , 025. 16/// ) 

COMPUTING  SATURATED  ACTIVITY  DISTANCE  TO  THE  FOURTH  AR 

AS*R**4 

D090  1=1,10 
R(I)=1.0*(I) 

Y(I)=Y(I)*R(I )**2.0 
90  CONTINUE 
H=1.0 

EXPODENTIAL  PORTION 

DO    200    1=10,80 
R(I)=1.0*(IJ 

Yd  )=C2*DEXP(-R(  I)/B)*R(  I)**2.0 

200    CONTINUE 

N  =  79 

SUM2=0.0D0 

SUM4=0.0D0 

M=N-1 

D0100    1=3, N, 2 

SUM2=SUM2+Y( I ) 
100    CONTINUE 

D0L10    1=2, M, 2 

SUM4=SUM4+Y< I ) 
110  CONTINUE 

S4=H/.3D1*(YU )+.4Dl*SUM4+.2Dl*SUM2+Y(N+l) ) 

WRITE(6,120)  S4 
120  F0RMAT(T19, « SIMPSONS  APPROXI  MATI  ON*  ,  D25  .  16///  ) 

AS(CD)*R**4 

H=1.0 
N=15 

DO  130  1=1,10 
R(I  )=1.0*< I) 
YCD(  I  J=YCD(I )*R(I)**2.0 
130  CONTINUE 

EXPODENTIAL  PORTION 

DO  260  1=10,80 
R(I)=1.0*( I) 
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c 
c 

c 
c 


YCD( I)=CD2*DEXP(-R(I )/BCD)*R(I )**2.0 
260  CONTINUE 

N=79 

SUM2=O.0D0 

SUM4=O.0D0 

M=N-1 

D0140  1=3, Nt 2 

SUM2=SUM2+YCD(I ) 
140  CONTINUE 

D0145  1=2, M, 2 

SUM4=SUM4+YCD( I ) 
145  CONTINUE 

SCD4=H/.3D1*(YCD< I ) +. 4D1*SUM4+ . 2D1 *SUM2+YCD( N+l ) ) 

WRITE(6,150)  SCD4 
150  F0RMAT(T19,» SIMPSONS  APPROX IMAT I  ON1 , 025 . 16/// ) 


COMPUTING  DIFFUSION  LENGTH 


AS2=S2 

AS4=S4 

ASCD2=SCD2 

ASCD4=SCD4 

MIG=AS4/(AS2*6.0I 

WRITE (6,290)  MIG 
2  90  FORM  AT ( T 19 , « M I GR  AT  I ON 

AGE=ASCD4/(ASCD2*6.0) 

WRITE(6,300)  AGE 
300  F0RMAT(T19,'AGE=»,D2  5.16///1 

DL=(MIG-( AGE+1.8) I**. 5 

WRITE(6,310)  DL 


310 


AREA=« , D25.16///I 


FORMAT (T 

STOP 

END 


19, 'DIFFUSION  LENGTH=«,D25.16///) 
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APPENDIX 


COMPUTER  PROGRAM  NO.  3 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  COMPUTE  AND  PLOT 
A  POLYNOMIAL  APPROXIMATION  OF  1/L**2  (K**2)   AS  A 

FUNCTION  OF  BORIC  ACID  CONCENTRATION  AND  TEMPERATURE. 
THE  BORIC  ACID  CONCENTRATION  IS  EXPRESSED  AS  THE 
MACROSCOPIC  ABSORPTION  CROSS-SECTION  OF  BORON.   THE 
TEMPERATURE  DEPENDENCE  FACTOR  (TN/TO)**.5  IS  BASED  ON 
A  REFERENCE  TEMPERATURE  OF  22.3  DEGREES  CENTIGRADE. 
REFERENCE  DATA  AT  THIS  TEMPERATURE  IS  FROM  REF.  (2). 


PROGAM  NOTATION 


EA  MACROSCOPIC  ABSORPTIN  CROSS-SECTION 

DL  DIFFUSION  LENGTH 

TN  MODERATOR  TEMPERATURE 

TO  MODERATOR  REFERENCE  TEMPERATURE 


REAL  DL,EA,B,C,U,V,SUM,TN,TO,A 

REAL  LABEL2/4H     / 

DIMENSION  DL(IO)  ,EA(10)  ,B(4) , SUM (4  I ,U(402)  ,V(401)  ,C(4, 

4),A(4J 

REAL*8  ITITLE  (12) 

RE AD (5,  10)  ((DL(I),EA(I)i,I=l,10) 
10  FORMAT(2F9.0) 

READ(5,34)  (  IT  I TL E (  I  )  ,  1  =  1 , 12  ) 
34  FORMAT (6A8) 

B(1)=0.0 

B(2)=0.0 

B(3)=0.0 

B(4)=0.0 

D025  1=1,10 

DL( I )=1.0/DL( I ) 

B(2)=B(2)+DL(I )*EA( I ) 

B(3)=B(3)+DL( I )*EA( I >**2 

B(4)=B<4)+DL( I )*EA( I ) **3 

25  CONTINUE 
D026  L=l,6 
SUM(L)=0.0 

26  CONTINUE 
D030  L  =  l,6 
D030  1=1,10 

SUM( L)=SUM(L)+EA(I )**(L) 
30  CONTINUE 

C(l,l)=10 

DO  40  M=l,4 

DO  40  N=l,4 

IF(M+N.EQ.2)  GO  TO  40 

C(M,N)=SUM(M+N-2) 
40  CONTINUE 

WRITE(6,50)  B( 1) , B ( 2 ) , B ( 3 » , B( 4 ) 
50  F0RMAK2X, 'VALUES  OF  B  «  ,  10X  ,  4F  1  5  .4  ) 

WRITE(6,60)  (SUM(L)  ,L  =  1,6) 
60  FORMAT (2X, •SUM  OFX I • , 1  OX , 6F 1 5 . 4 ) 

WR I TE ( 6 , 7  0 )  ( i  C (  M , N )  , M= 1 , 4 )  , N  =  1 , 4 ) 
70  F0RMAT(2X, 'VALUES  OF  C',16F12.4) 
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80 


90 


91 


92 


93 


94 


CALL 

WRIT 

FORM 

U(l) 

DO  9 

V(  I) 

U(  1  + 

CONT 

CALL 

,LAS 

CALL 

LAST 

TN=3 

T0=2 

A(l) 

A(2) 

A(3) 

A(4) 

A(4) 

A(3) 

A(2) 

A(l) 

DO  9 

V(  I) 

U(  i* 

CONT 

CALL 

LAST 

TN  =  3 

A(l) 

A(2) 

A(3) 

A(4) 

A(4) 

A(3) 

A(2) 

A(l) 

DO  9 

V(  I  ) 

U(  1  + 

CONT 

CALL 

LAST 

TN=3 

A(  1) 

A(2) 

A(3) 

A(4) 

A(l) 

A(2) 

A(3) 

A(4) 

DO  9 

V(I) 

U(  1  + 

CONT 

CALL 

LAST 

TN  =  3 

A(  1) 

A(2J 

A(3) 

A(4) 

A(l) 

A(2) 

A(3) 

A(4) 

DO  9 

V(  I) 

U(  1  + 

CONT 

CALL 


B,4, 

B(1J 

,OLUT 


K  S  ) 

,B(2) ,B(3) ,B(4) 

ION' ,6X,4F10.4) 


1 

I)*(B(2)+U(I)*B(3J+(U(I )**2)*B(4) ) 

.001 


t  EA, 
ltU, 


DL, 1, 3, LAB EL2, I  TITLE, 0,0,0, 0,0, 0,6, 9,1 
V,2,0,LABEL2,ITITLF,0,0,0,0,0,0,6,9,1, 


SIMQ(C, 
E(6,80) 
AT(2X,  'S 
=  .01 

0  1=1,40 
=  B  (  1 )  4  U  ( 

1)=U(I  )4 

INUE 

DRAW( 10 
TJ 

DRAW(40 
) 

08.0 
96.15 
=  B(1) 
=  B(2) 
=  B(3) 
=  B(4) 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 

1  1=1,40 
=A( 1)+U( 
1)=U( I )4 
INUE 

DRAW  (401  ,U,V,2,0,LABEL2,ITITLE,0,0,0,0,0,0,6,9,1, 
} 


1 
I)*< 

.001 


5)*A(4)*.8862 
5)*A(3)*.8862 
5)*A(2)*.8862 
5)*A(1)*.8862 

A(  2)4U(I  )*A(3)+(U(  I  )**2)*A(41  ) 


28.0 
=  B(1) 
=  B(2) 
=  B(3) 

=  B(4) 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 

2  1=1,40 
=A( 1)+U( 
1)=U( 1 )4 
INUE 

ORAW(40 
) 

38.0 
=  B(1) 
=  B(2) 
=  B(3) 
=  B<4) 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 

3  1=1,40 
=A(1)+U( 
1)=U(I  )4 
INUE 

DRAW (40 
) 

43.0 
=  B(1) 
=  B(2J 
=  B(3) 
=  B(4) 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 
=( (TO/TN 

4  1=1,40 
=  A(1H-U( 
1)=U( I  )4 
INUE 

DRAW (40 


)** 

1 

I)*(A(2)+U(I)*A(3)+(U<I )**2)*A(4) ) 

.001 


5)*A(4)*.8862 
5)*A(3)**8862 
5)*A(2)*.8862 
5)*A( 1  J*. 8862 


ltUt V,2t Of LABEL2, 1  TITLE t 0,0,0, 0,0, 0,6, 9,1, 


)**.5)*A( 1 )*.8862 

)**.5)*A<2)*.8862 

)**.5)*A(3)*.8862 

)**.5)*A(4)*.8862 

1 

I )*( A(21+U(I )*A(3)+(U(I)**2)*A(4) ) 

.001 

1,U,V,2,0,LABEL2,ITITLE,0,0,0,0,0,0,6,9,1, 


)**. 
)**. 

)**. 


5)*A( 1 )*.8862 
5)*A(2)*.8862 
5)*A(3)*.8862 

5)*A(4)*.8862 


1 

I)*(A(2)+U(I)*A(3)+(U(I  )**2)*A(4)  ) 

.001 

1,U,V,3,0,LABEL2,ITITLE,0,0,0,0,0,0,6,9, 1, 
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LAST) 

STOP 

END 
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